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ABSTRACT 

Past and future dynamical evolution of all 64 long period comets having l/a or ;< 1 x 10~ 4 au~' 
and g osc > 3.0 au and discovered after 1970 is studied. For all of them we obtained a new, ho- 
mogeneous set of osculating orbits, including 15 orbits with detected non-gravitational para- 
meters. The non-gravitational effects for eleven of these 15 comets have been determined for 
the first time. This means that more than 50% of all comets with perihelion distances between 
3 and 4 au and discovered after 1970 show detectable deviations from a purely gravitational 
motion. Each comet was then replaced with the swarm of 5001 virtual comets well represent- 
ing observations. These swarms were propagated numerically back and forth up to the 250 au 
heliocentric distance, constituting sets of original and future orbits together with their uncer- 
tainties. This allowed us to show that the 1 /a OI i distribution is significantly different in shape 
as well as the maximum position when including non-gravitational orbits. Next we followed 
the dynamical evolution under the Galactic tides for one orbital revolution to the past and 
future, obtaining orbital elements at previous/next perihelion passages. We obtained a clear 
dependence of the last revolution change in perihelion distance with respect to the 1 /a 0I i, what 
confirmed theoretical expectations. 

Basing on these results we discuss the possibility of discriminating between dynamically 
new and old comets with the aid of their previous perihelion distance. We show that about 
50 per cent of all investigated comets have their previous perihelion distance below the 15 au 
limit. This resulted in classifying 3 1 comets as dynamically new, 26 as dynamically old and 
7 having unclear status. We showed that this classification seems to be immune against to 
perturbations from all known stars. However discoveries of new, strong stellar perturbers, 
while rather improbable, may change the situation. 

We also present several examples of cometary motion through the Jupiter-Satum bar- 
rier, some of them with the previous perihelion distance smaller than the observed one. New 
interpretations of the long period co mets source pathways are also discussed in the light of 
suggestions of iKaib & Ouinnl J2009t) . 
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1 INTRODUCTION 

As Martin Duncan d2009h recently pointed out, the forthcom- 
ing deep, wide-field observational surveys (both ground-based and 
space-based) soon be providing completely new amount of data on 
very large perihelion distance comets, far behind Saturn. These may 
extend our possibilities to investigate the source regions of the long 
period comets (LPCs). But up to the date, the only observational 
constrains on the LPCs source is the detailed analysis of their ob- 
servations, ranging up to 10 au or so, determining their orbits and 
studying thei r past dynamical evolution. Contin uing our effort in 
this field (see lKrolikowska & Dyb czynski 2oTol. hereafter Paper I) 
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we analysed the sample of next 64 comets from the so called Oort 
spike, all of them having the osculating perihelion distance greater 
than 3 au, wherein four comets with q osc > 3.0 au and determinable 
non-gravitational (here and after NG) orbits were taken from Paper 
I. The detailed description of our sample can be found in Section 
[2] In short, the reason for such a selection was to deal with comets 
with negligible NG forces due to their large heliocentric distance. 
However, during the data analysis we attempted to determine the 
NG forces parameters and it appeared that for another 1 1 comets it 
was reasonable to use such a model. Thus, including four comets 
from Paper I we have in total 15 large perihelion distance comets 
with NG parameters obtained and used in the original/future orbit 
determinations, details can be found in Sections |2. 1 1 — 12.31 In this 
paper we again ask for the source region of the Oort spike comets. 
To this purpose we carefully analysed past and future dynamical 
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evolution of 64 comets under planetary and Galactic perturbations, 
see Section [3] This, additionally, makes an opportunity to observe 
how the mechanism widely called Jupiter-Saturn barrier works in 
practice, see Sections Q7T] and [3] The widely disputed problem of 
discriminating between dynamically (and physically) new and old 
comets is revisited in Section [4] including a discussion of our res- 
ults in th e light of the newly p roposed alternative cometary origin 
scenario dKaib & Quin r] |2009h . Finall discussion and conclusions 
are presented in Section|5] 



1.1 On the Jupiter-Saturn barrier 

The concept of Jupi ter- Saturn b arrier can be traced back by 30 
years or so. In 1981 IFernandezI presented the dependence of the 
comet energy changes due to the planetary perturbation on the peri- 
helion distance. It appeared, that for perihelion distances smaller 
than 15 au this perturbation is comparable with LP comet binding 
energy. 



Later on, as a result of some numerical simulations IWeissmanl 

d 19851) stated that 65 per cent of comets coming closer to the Sun 
than Saturn (and 94 per cent for Jupiter) will be ejected from the 
Solar System as a result of planetary perturbations. 

These results of planetary perturbations investigations 
were then combined wit h the model of Galactic disk tides. 
iMatese & Whitman! dl989t> showed that it is possible to calculate 
the minimum cometary semimajor axis, for which Galactic per- 
turbations can decrease cometary perihelion distance from above 
the strong planetary perturbation border down below the observ- 
ability limit in one orbital period. Using 15 au as the former and 
5 au as the later they obtained the minimum semimajor axis to 
be equal 20 000 au. Inves tigating Galactic perturbations in a way 
somewhat similar to ours, lYabushitald 19891) obtained a very strong 
dependence of the perihelion distance reduction by Galactic tides 
on the cometary semimajor axis, namely Aq ~ a 6 - 3±0 - 2 . His con- 
clusion was that 25 000 au is the most probable semimajor axis of 
LP comets arriving at the vicinity of t he Su n for the first time. 
Finally, in a review paper, IFernandezI d 19941) devoted a separate 
section to the description of the Jupiter-Saturn barrier mechan- 
ism. This w a s also d iscussed by m a ny au th ors, see for e xampl e 
FestouetalJ Jl993h: iLevison et al.l feOOll) : bones et al.l d2004l) : 



Fernandez! d2005l) ; lMorbidellil d2005l) . 



This effect can be described in short as follows: the Galactic 
perturbations can result in a continuous cometary perihelion drift 
towards the Sun. The rate of this drift strongly depends on cometary 
semimajor axis, which is almost constant here. If we agree that 
the perihelion distance smaller than 10-15 au results in the ejection 
from the Solar System by perturbations from Jupiter and Saturn the 
only possibility for a LP comets from the Oort cloud to become 
observable is to have their perihelion distance reduced from above 
10-15 au down below the observability limit, say 3-5 au. Basing 
on the Galactic tide model it is possible to calculate the minimum 
necessary semimajor axis to accomplish this. 

Such a calculation can be performed in several slightly differ- 
ent formalisms but only two factors can change the result signific- 
antly: the amplitude of the necessary perihelion distance reduction 
and the assumed Galactic disk matter density p . In almost all papers 
the authors demand the perihelion reduction by -10 au (from above 
10-15 au to below 3 -5 au) but they use differen t disk matter densit- 
ies. In earlier papers dMatese & Whitman! 1989llFesiou et al.ll993l) 
they use d p = 0.185 solar masses per cubic parsec, la ter (see for 
example: Iwiegert & Tremainelll999l: lMorbidellill2005l) p = 0.150 
was used. In all recent papers the value of p = 0.100 solar masses 



per cubic parsec is used, see lLevison et al.l d200ll) for a supporting 
discussion of this density value. 

As a consequence, while previously obtained results (a >~ 
20000 au) were roughly in line with o bservations, now the situation 
is more complicated. In a review by iDones et al.l d2004l) one can 
read: 'If we assume that a comet must come within 3 au of the Sun 
to become active and thus observable, Aq must be at least ~10 au - 
3 au = 7 au. It can be shown that, because of the steep dependence 
of Aq on a, this condition implies that a > 28000 au'. But there are 
many observed Oort spike comets with much smaller semimajor 
axes! 

There is also an additional aspect of the Jupiter- Saturn bar- 
rier. The resulting semimajor axis threshold value is often used as 
a definition of the border between the outer (observable) and in- 
ner (unobservable) parts of the Oort cloud. With the current value 
(a > 28000 au ) this seems to b e prob lematic. 

Recently Kaib & Ouind 12009?) revisited the Jupiter- Saturn 
barrier problem and give several interesting new conclusions. We 
will discuss them in Section l4"l4l 

In what follows we present some examples showing how dif- 
ferent from the above theory is the behaviour of many observed 
Oort spike comets. 



2 THE SAMPLE OF LARGE PERIHELION LP COMETS 

The latest published Catalogue of Cometary Orbits 
dMarsden & Williamsll2008l hereafter MWC08) includes 154 Oort 
spike and hyperbolic comets (l/a ri < 10~ 4 au~') with orbits 
determined with the highest precision, i.e. with the quality class 1 , 
according to the classification introduced bv lMarsden et al.l dl978l ). 
Their 1 /a or \ distribution, based on catalogue data, is shown in 
FigQJa). In the present investigation we chose to study a sample of 
the Oort spike comets with large perihelion distances, q ^ 3.0 au. 
The MWC08 contains 76 such comets. According to the availabil- 
ity of astrometric data we restrict this sample to comets discovered 
after 1970. Additionally, to have orbits definitely determined we 
did not analyse three still potentially observable comets in October 
2010: C/2005 L3, C/2006 S3 and C/2007 Dl. Thus, our sample 
was reduced to 62 comets which orbits are of quality class 1 
in MWC08. Recently, this sample has increased by two comets 
(C/2006 YC, C/2007 Yl) that have orbits with the quality class 2 
in MWC08 but now their number and interval of observations 
have increased significantly allowing for much more precise orbit 
determination. This makes the total number of 64 LP comets 
studied in the present paper. The combined l/a or ; distribution 
containing MWC08 comets superposed and supplemented by our 
results is presented in Fig. [T](b). One can observe the influence 
of our new results on the overall shape of the histogram as well 
as position of l/a ri -maximum. The black part consists of all 
86 comets investigated in Paper I and present research. It should be 
stressed here, that during the orbit determination process we obtain 
l/ a ori value as well as its uncertainty. This uncertainties are taken 
into account when constructing the black part of the histogram 
depicted in Fig.Q](b). It is evident, that incorporating NG effects in 
orbit determination (where possible) moves the overall distribution 
towards smaller semimajor axes. 

For each comet from the sample we determined its osculat- 
ing nominal orbit (gravitational and NG if it was possible; see next 
section) based on the data selected and weighted according to the 
methods described in great detail in Paper I (see also section 12.11 
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Figure 1. Distribution of \/a ml x 10 6 au -1 for all LP comets with class 1 orbits and l/a or j < 10~ 4 au~' . (a): 153 comets with l/a or j taken directly from the 
MWC08 (splitting product C/1996 Jl-A omitted) , (b): 165 Oort spike comets; black part of the histogram denotes 86 comets with orbits determined in this 
paper and Paper I, We increased the MWC08 sample by 12 comets and determined the NG orbits for 37 comets. 



Table 1. Original and future semimajor axes derived from pure gravitational nominal solutions (columns 3^) and NG nominal solutions (columns 5-6) for 15 
large perihelion distance comets with detectable NG effects; the number of NG parameters determined for NG solutions is given in the column 1 1 . The rms's 
and number of residuals are given in the columns 7-8 and 9-10, respectively. The numbering in column 1 is consistent with the chronology of the discovery 
of all 64 considered comets. Solutions for comets C/1997 BA 6 , C/1997 J2, C/1999 Yl and C/2000 SV 74 are taken from Paper I. 
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10 
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C/1974 Fl Lovas 


+36.82± 


2.50 


+518.11± 


2.50 


+40.48± 


5.00 


+522.07± 6.81 


1.09 


1.08 


273 


273 


2 


10 


C/1980 El Bowell 


+27.52± 


1.24 


-16014.± 


1. 


+53.35± 


3.87 


-16011.± 


4. 


1.17 


1.06 


388 


387 


2 


11 


C/1983 01 Cernis 


+47.79± 


1.58 


-191.39± 


1.58 


+60.83± 36.12 


-186.95± 


2.08 


1.14 


1.11 


461 


461 


2 


12 


C/1984 W2 Hartley 


+13.89± 


22.32 


-27.92± 


9.62 


+20.25± 


8.56 


-31.55± 


8.57 


1.89 


1.87 


107 


107 


1 


21 


C/1997 BA 6 Spacewatch 


+1.49± 


0.35 


+371.77± 


0.35 


+31.83± 


1.15 


+402.48± 


1.72 


0.74 


0.67 


1054 


1054 


3 


22 


C/1997 J2 Meunier-Dupouy 


+38.83± 


0.57 


-2.72± 


0.57 


+44.64± 


0.88 


+14.72± 


0.91 


0.67 


0.53 


2881 


2863 


2 


25 


C/1999 H3 LINEAR 


+65.99± 


0.71 


-8.89± 


0.71 


+124.66± 


3.88 


-10.03± 


1.05 


0.70 


0.51 


1739 


1722 


3 


32 


C/1999 Yl LINEAR 


+42.92± 


0.88 


+350.42± 


0.88 


+47.35± 


0.94 


+345.69± 


1.46 


0.61 


0.48 


1747 


1749 


3 


34 


C/2000 CT 54 LINEAR 


+38.75± 


1.35 


+557.96± 


1.35 


+72.86± 


3.06 


+585.28± 


2.53 


0.90 


0.75 


418 


417 


2 


37 


C/2000 SV74 LINEAR 


+50.23± 


0.40 


-85.55± 


0.40 


+92.31± 


0.85 


-54.88± 


0.60 


1.11 


0.71 


4389 


4349 


3 


47 


C/2002 R3 LONEOS 


+38.98± 


0.57 


+9.88± 


0.57 


+48.16± 


3.07 


+3.86± 


6.33 


0.55 


0.52 


2533 


2530 


3 


54 


C/2005 Bl Christensen 


+7.14± 


0.33 


+234.38± 


0.33 


+3.99± 


0.61 


+239.09± 


0.47 


0.45 


0.43 


2991 


2985 


3 


55 


C/2005 EL173 LONEOS 


+45.70± 


0.47 


-35.29± 


0.47 


+44.81± 


0.99 


-19.12± 


0.91 


0.47 


0.36 


631 


632 


2 


57 


C/2005 Kl Skiff 


+7.66± 


1.22 


-79.34± 


1.22 


+11.01± 


2.75 


-82.36± 


3.35 


0.52 


0.52 


1257 


1254 


2 


61 


C/2006 S2 LINEAR 


+73.92± 


2.75 


-30.83± 


2.75 


+72.52± 


8.14 


-10.77±1 


8.00 


0.45 


0.44 


346 


346 


2 



therein). This allows us to construct a homogeneous sample of 
cometary osculating orbits which are the starting point for us to ob- 
tain the original and future orbits with their uncertainties and then 
to study cometary dynamical evolution under the Galactic tides. 

Our set of large perihelion LP comet orbits with l/a or i < 
10~ 4 au~' contains only one comet with slightly negative value 
of l/fl ori , namely C/1978 G2. The MWC08 attributes original 
hyperbolic orbits to three more comets: C/1942 C2, C/2002 R3, 
C/2005 Bl. The first was discovered before 1970, therefore it was 
excluded from our sample. In the present paper, orbits of two other 
comets are determined based on a significantly longer interval of 
observations than in MWC08 and now their original (gravitational) 
orbits derived by us are elliptical. Later in this paper, it appears that 
comet C/1978 G2 may also be of local origin. Therefore, we call 
this sample the large perihelion distance Oort spike comets. For 
comets with small perihelion distances we have an essentially dif- 
ferent situation. In MWC08 nearly 20 comets with q < 3.0 au have 



negative values of l/a ri ( see FigQJ). It turns out that for these 
comets important is the role of NG effects in the process of de- 
termining their osculating orbits: most of small perihelion comets 
on hyperbolic original orbits in the gravitational case have elliptical 
original orbits when the NG forces are taken into account (see next 
section). 

Observed perihelion distance and ecliptic inclination distribu- 
tions of all investigated comets are also shown and discussed in 
Sectionl4~T1 

2.1 Detection of the non-gravitational forces in large 
perihelion distance LPCs 

It is well-known that the actual motion of the comet is not purely 
gravitational (here and after GR). This means, that very small 
uncertainties of orbital elements in a GR model are often only 
the formal errors calculated in the process of orbit determination. 
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Table 2. The past distributions of swarms of VCs in terms of returning [R], escaping [E], including hyperbolic [H] VC numbers for dynamically old comets. 
Aphelion and perihelion distances are described either by a mean value for the normal distributions, or three deciles at 10, 50 (i.e. median), and 90 per cent. 
In the case of mixed swarm the mean values or deciles of Q and q are given for the returning part of the VCs swarm, where the es cape limit of 1 20 000 au was 
generally used. The upper-a index in columns 5-7 means that this part of mixed swarm includes the nominal orbit. For comparison we included the osculating 
perihelion distance in the third column; in the fourth column - the galactic latitude of the perihelion direction is given. The last two columns present the value 
of 1/flori and the percentage of VCs, that we can call dynamically new, based on previous q statistics. Fifteen comets with NG effects including four objects 
from the Paper I (No: 21,22,32,37) are indicated by upper-NG index located behind the comet designation (column 2). 
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1 


C/1972 LI 


4 28 


40.3 


5001 








33.9 - 39.0 - 46.5 


4.06-4.18-5.08 


51. 1± 


6.2 


n s - n 9 - ft 1 


2 


C/1973 Wl 


3.84 


-48.7 


5001 








22.7 - 27.5 - 35.2 


4.14-4.47-5.53 


72.6±12.3 


0.7 - 0.4 - 0.2 


3 


C/1974 F1 NG 


3 01 


-21.7 


5001 








42.7 - 49.4 - 58.6 


6.18-8.98- 16.46 


40.5± 


5.0 


J7.4 - 1J.J - J.I 


5 


C/1976 D2 


6 88 


44.5 


5001 








30.2- 35.1 -42.2 


4.76 - 5.50 - 6.00 


56.9± 


7.3 


0-0-0 


6 


C/1976U1 


5 86 


30.5 


4556° 


445 


81 


26.7 - 41.8 - 78.4[R] 5.51 - 5.72 - 56.21[R] 


45.6±21.9 


19 9 - 97 1 - 94 1 

JZ.Z. iit.J 


9 


C/1979 M3 


4.69 


-14.2 4781" 


220 


13 


33.2 - 47.9 - 79.1[R] 4.30 - 4.48 - 24.35[R] 


41.1±14.5 


9ft Q - 1 ft Q - 1 S ft 
z.u .y iu.7 u.v/ 


10 


C/1980 E1 NG 


3.36 


13.3 


5001 


() 





34.3 -37.5-41.3 


1.82-2.16-2.41 


53.3± 


3.9 


0-0-0 


11 


C/1983 01 NG 


3.32 


-55.2 4503" 


498 


210 


18.2-30.9-68.6 


3.53 - 4.74 - 52.74 


60.8±36.1 


32.8 -27.6 - 24.8 


13 


C/1987 Fl 


3.62 


26.5 


5001 








31.0-34.4-38.5 


4.82 - 5.40 - 6.39 


58.2± 


5.0 


0.1-0-0 


14 


C/1987 HI 


5.46 


17.6 


5001 








40.8 - 44.0 - 47.8 


8.66-9.79- 11.53 


45.5± 


2.8 


42.5 - 0.3 - 


IS 


C/1993 Fl 


5.90 


45.1 


5001 








28.3-31.9-36.7 


7.23 - 8.00 - 9.48 


62.6± 


6.3 


5.8-0.1-0 


22 


C/1997 J2 NG 


3.05 


0.5 


5001 








44.8 ± 0.9 


2.801 ± 0.017 


44.6± 


0.9 


0-0-0 


24 


C/1999 F2 


4.72 


46.0 


5001 








39.1 -42.6-46.8 


6.21 - 7.03 - 8.59 


46.9± 


3.3 


2.3-0.1-0 


25 


C/1999 H3 NG 


3.50 


49.9 


5001 








16.0 ± 0.5 


3.629 ± 0.014 


124.7± 


3.9 


0-0-0 


28 


C/1999 N4 


5.51 


23.1 


5001 








29.4 ± 0.7 


6.527 ± 0.092 


68.2± 


1.7 


0-0-0 


29 


C/1999 S2 


6.47 


-22.3 


5001 








32.6-35.4-38.8 


8.17-8.79-9.77 


56.4± 


3.8 


6.6-0-0 


30 


C/1999 Ul 


4.14 


-25.0 


5001 








48.0-52.8-58.8 


6.23 -7.84- 11.39 


37.9± 


3.0 


18.6-2.5-0 


32 


C/1999 Y1 NG 


3.09 


-62.7 


5001 








42.2 ± 0.8 


5.82-6.11 -6.46 


47.4± 


0.9 


0-0-0 


34 


C/2000 CT^ 4 G 


3.16 


-36.7 


5001 








27.4 ± 1.2 


3.88 - 4.04 - 4.24 


72.9± 


3.1 


0-0-0 


36 


C/2000O1 


5.92 


-23.5 


5001 








34.4 - 38.5 - 43.5 


7.16-7.87-9.19 


52.0± 


4.7 


4.1-0-0 


37 


C/2000 SV^f 


3.54 


21.7 


5001 








21.7 ±0.2 


3.791 ± 0.009 


92.3± 


0.8 


0-0-0 


45 


C/2002 J5 


5.73 


33.0 


5001 








33.77 ± 0.39 


7.966 ± 0.098 


59.2± 


0.7 


0-0-0 


47 


C/2002 R3 NG 


3.87 


-45.8 


5001 








38.3-41.5-45.1 


4.69-5.15-5.94 


48.2± 


3.1 


0-0-0 


55 


C/2005 EL^73 


3.89 


-23.2 


5001 








43.66 ± 0.99 


7.82-8.31 -8.89 


44.8± 


1.0 


0-0-0 


61 


C/2006 S2 NG 


3.16 


-14.9 


5001 








24.1 -27.5 -32.2 


3.193 -3.221 - 3.294 


72.5± 


S.l 


0-0-0 


63 


C/2007 JA 2 i 


5.37 


31.4 


5001 








30.75 ± 0.95 


5.86 - 5.94 - 6.04 


65. 1± 


2.0 


0-0-0 



These small uncertainties, in turn, may lead us to believe that we 
know the true orbit of a comet so perfectly. When incorporating 
NG effects in the dynamical model we often obtain a slightly larger 
formal uncertainties but we believe that the resulting orbit better 
represents the real cometary motion, especially when extrapolating 
out of the observed time interval. In particular, this may be the case 
of some near-parabolic orbit when taking into account NG effects 
may even lead to changes in the sh ape of the original baryce ntric 
orbit from hyperbolic to elliptical d Krolikow ska 200 1[ I2OO6I and 
Fig. 5 in Paper I). It has been shown that almost all comets with 
hyperbolic original GR orbits have original NG orbits elliptical. 

In Paper I, our sample of comets with the NG effects was 
composed of objects with NG orbits characterized by a clear de- 
crease in the rms compared to their rms for purely gravitational 
orbits. We found then 26 comets with NG orbits very well de- 
termined for which original reciprocals of semimajor axes in the 
NG model, l/fl ri.NG> were located within the range from zero to 
100 x 10~ 6 au~ 1 . In preparing that sample many comets of a purely 
GR orbit in Oort spike has been excluded because their NG orbit 
gave 1 /a or i ng > 100 x 10~ 6 au -1 . That sample of 26 comets con- 
tained only four objects with perihelion distances exceeding 3 au. 

In the present study we approach the issue differently. Em- 
phasis is laid on the completeness of the sample. This allowed 
us not only to investigate the dynamic behaviour in the previous 



and future perihelion but also to make a statistical analysis of the 
observed sample of Oort spike comets and their apparent source 
region. Knowing that the NG effects in the motion of LP comets 
are hardly determined we concentrated on a sample of comets with 
large perihelion distances in order to keep the highest possible orbit 
reliability even for objects with indeterminable NG effects. 

However, we have not given up on determining the NG effects 
in these comets though, for comets with large perihelion distances 
the NG effects are determined even more difficult than for comets 
approaching closer to the Sun. Therefore, we decided to set out 
the NG orbit regardless of whether it is a significant drop in rms. 
We assumed that even if this decrease is negligible, the model of 
the NG motion should give a better representation of the actual 
comet's orbit and its uncertainty (NG model consists from one to 
four additional parameters, so the uncertainties of orbital elements 
can be significantly larger). Once determined, the NG parameters 
with reasonable accuracy, we recognized the model as more real- 
istic than the purely GR model. The obtained NG models were also 
reviewed for O-C time variations and O-C distributions (a more de- 
tailed description of the methods used are given in Paper I). In this 
way, from the sample of 64 comets, we have determined NG orbits 
for 15 objects with perihelion distances greater than 3.0 au. It ap- 
peared that all of them have q lower or equal to 4 au. Comparison 
between GR and NG solutions for these fifteen comets is given in 
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TableQ]where comets are ordered by discovery date. Four comets 
from Paper I, with a very well-determined NG effects are included 
here. Finally, in our sample, we have 23 objects with perihelion dis- 
tance in the range of 3.0 au < q < 4.0 au (see also Table[7J. It means 
that only eight comets of this perihelion distance range still seem 
to have undetectable NG effects. It should be summarized that it 
is clearly more than a half of all comets with perihelion distances 
between 3 and 4 au (discovered after 1970) that show small devi- 
ations from purely gravitational motion, which can be detectable 
either through a decrease in rms for NG model of motion or by 
improvements when analysing the differences in OC distribution 
and/or OC time variations between GR and NG models. 

To determine the NG cometary orbit we used t he same form- 
alism as in Paper I that was originally proposed bv lMarsden et al.l 
dl973l) . This formalism introduces three orbital components of the 
NG acceleration {A\,Ao, and A3, i.e. radial, transverse and normal 
components, respectively) acting on a comet in a case of sublima- 
tion symmetric relative to perihelion. The asymmetric NG model 
introduces additional parameter T, the time displacement of the 
maximum of the g(r(/ — t)). From orbital calculations, the NG para- 
meters Ai,Ax, and A3 and eventually z should be derived together 
with six keplerian orbital elements with in a given observat ional 
time interval (more details are given in lKrolikowskall2006l) . We 
re alize that the standard g (r)-function used by us has been obtained 
bv lMarsden et al. U 19731) phenomenologically on the basis of wa- 
ter sublimation. We have tried to determine the NG orbits for the 
investig ated comets assumin g more general form of g(r)-like func- 
tion (see Krolikowska 2004), however, it appeared that the rms and 
OC time variations were practically the same. So we came to the 
conclusion that it is best to use standard form of the g(r)-function. 
Past experience in determining the NG orbit on the basis of various 
g(r)-like functions jKr61ikowskall2004l) justifies such an approach. 

For six comets studied here we were able to derive all three 
parameters of NG acceleration, for the next eight - the radial and 
the transverse component of the NG acceleration, in one case 
(C/1984 W2) - only the radial term. 



2.2 Original and future orbits. Creating swarms of virtual 
comets. 

The very first step in studying past and future motion of LP comets 
is to determine their original and future orbits. As usual, we call an 
orbit original when traced back out of reach of planetary perturb- 
ations (assumed in this paper to happen at 250 au from the Sun). 
Similarly, the future orbit can be obtained after following the mo- 
tion of a comet up in time as far as it reaches the same heliocentric 
distance of 250 au. 

To derive original and future orbital parameters (including se- 
mimajor axes) as well as their uncertainties we have examined the 
evolution of thousands of virtual comet (VC ) orbits using t he Sit- 
arski's method of the random orbit selection (Sitarski 1998). More 
details are given in Paper I. With this method, we construct oscu- 
lating swarms of comets that follow the normal distribution in the 
space of orbital elements and eventually NG parameters (the 6- 
9 dimensional normal statistics). Similarly as in Paper I we fill a 
confidence region with 5 000 VCs for each nominal solution (GR 
and NG if determined) for each considered comet. 

Next, each of VC from the osculating swarm was followed 
numerically from its position at osculation epoch backwards and 
forwards until this individual VC reached a distance of 250 au from 
the Sun. The equations of comet's motion have been integrated nu- 
merically using the recurrent power series method dSitarskilll98^. 
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Figure 2. Shifts of l/o ri due to the NG effects for 15 of investigated comets 
and all comets from Paper I. 



l2002h . taking into account perturbations by all the planets and in- 
cluding the relativistic effects. In this way we were able to obtain 
the nominal original and future barycentric orbits of each comet as 
well as the uncertainties of the derived values of orbital elements by 
fitting t he normal distributi on to each original and future cometary 
swarm d Kro likowska 2001). Obtained original and future reciproc- 
als of semimajor axes with their uncertainties are given in Tables|2] 
-|4]and Tables 1914101 where comets are ordered by discovery date. 
As it was mentioned before, only one comet in the investigated 
sample, C/1978 G2, formally have a negative l/a OI i but large un- 
certainty of this value clearly does not exclude Solar System origin 
of this comet. 

The original and future swarms of VCs became the starting 
data to study the dynamical evolution of each individual comet un- 
der the Galactic tides (Section|3}. 



2.3 Differences between non-gravitational and gravitational 
models 

For 15 comets with the determined NG orbits we present only the 
results for the dynamic evolution of NG swarms of VCs (Tables f2] 
-|4]and Tables [9l410t. However, we derived also their original and 
future orbits starting from the GR osculating orbits. Differences 
between the results obtained we will discuss in terms of differ- 
ences in values of l/fl ri an d V fl fut> an d possible differences in 
the assessment of the dynamical status of investigated comets (Sec- 
tion|4j- In this section we will focus on the first issue. 

Figure|2]shows the differences between original reciprocals of 
semimajor axes for GR and NG models. Comets investigated here 
are presented by filled symbols (15 comets with q ^ 3.0au), while 
the open symbols represent 22 comets with q < 3.0 au from Paper I. 
When drawing conclusions, remember that comet sample analysed 
in the Paper I was constructed on the other principles than a sample 
of comets currently being examined. However, Fig.|2]clearly shows 
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Figure 3. Previous perihelion distance - eccentricity distribution for 
C/2000 SV74. All 5001 VCs were stopped at the previous perihelion. The 
centre of the dashed circle marks the nominal orbit. 
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Figure 4. Previous perihelion distance - eccentricity distribution for 
C/2005 Ql. All 5001 VCs were stopped at the previous perihelion. The 
centre of the dashed circle marks the nominal orbit. 



that in the case of comets with low q (less than, say, q ^ 2 au), the 
differences A(l/a or i) = 1 /a or ;,NG — l/ fl ori.GR Me much greater than 
that for comets with large q {q greater than 3 au). For more than 
50 per cent of large perihelion comets A(l/a or i) < 20 x 10~ 6 au 
_1 . It is worth noting that of 15 comets with l/a or i.GR within 
Oort spike, only one has an NG orbit slightly outside the peak 
(C/1999 H3, l/flori,NG = (124.66±3.88) x 10~ 6 au " 1 ) while hav- 
ing the greatest value of A( 1 /a or ; ) = 58 x 10~ 6 au ~ 1 resulting from 
planetary perturbations. When constructing the sample of NG Oort 
spike comets in Paper I, we obtained several such cases. 

Usually the uncertainties of orbital elements including original 
and future semimajor axes are larger for NG orbits than GR or- 
bits. Thus, initial NG swarms of VCs for original and future orbit 
calculations are more dispersed than GR swarms. In our opinion, 
the NG swarms of VCs better reflect actual knowledge of studied 
cometary orbits. This affects the determination of 1 /a ri,NG uncer- 
tainty, which is generally a factor of 2-3 larger than the uncertainty 
of l/a ri,GR determination for the same comet, as it can be ob- 
served in TableQ] An extreme case is 1 /a ri,NG = (60.8 ± 36. 1) x 
10" 6 au" 1 for comet C/1983 01 while the purely GR swarm gives 
l/ a ori,GR = (47.8 ±1.6) x 10~ 6 au -1 , which is solution formally 
more than one order of magnitude more accurate. 

For the remaining 49 comets NG effects are indeterminable. 
Hence, we present the results for the purely gravitational swarms 
of VCs. However, only 8 of these comets have perihelion distances 
less than 4 au. With the apparent trend of significant reduction of 
differences between NG orbit and GR orbit with the increasing 
perihelion distance (Fig. [2} it seems reasonable that the vast ma- 
jority of real orbits of these 49 comets is well-targeted despite of 
omitting the NG effects. 



3 PREVIOUS AND NEXT PERIHELION PASSAGES 

Starting from original or future LPCs orbits we followed their 
dynamical evolution under the Galactic tides, where both disk 
and central terms were included. This i s possible on l y in t he ab- 
sence of any other perturbing forces. iDvbczvhskil J2OO6I) have 
shown, that none of the known star influenced the motion of LP 
comets significantly in the last, say 10 million years, and the 
same holds for the same interval in the future. The similar con- 
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Figure 5. Past and future heliocentric distance changes of all VCs repres- 
enting C/2002 J5. Grey (cyan) lines present heliocentric distance changes 
of all VCs except the nominal ones, for which black line is used. Addition- 
ally the distribution of the future VCs swarm (black histogram) in time is 
presented. 



elusion one can find in Delsemme Jl987h. iMatese & Whitman! 
1992I). IWiegert & Tremain 3~ifl999h. Matese & Lissauen 1 d2004h . 



Dybcz ynski ( 2 002b, Emel'vanenko et all d2007l) or recently 
Kaib & Ouinnl d2009T ). This time interval is comparable with the 
orbital period of a comet having the semimajor axis of 50000 au 
therefore we decided to follow the motion of each comet for one 
orbital period to the past and future. Additionally, longer numerical 
integrations would show rather artificial cometary motion in cases 
of previous/next perihelion in planetary region since the planet- 
ary perturbations cannot be taken into account for obvious reasons. 
More discussion on the influence of omitting stellar perturbations 
on our results one can found in Section [5] In the way described 
above, we obtain previous and next perihelion passage distances or 
detect past or future escape, see Table [5] The final orbits obtained 
from this calculations we call previous and next . 

In this paper we call a comet (or more precisely each indi- 
vidual VC) as returning [R] if it goes no further than 120000 au 
from the Sun. All other comets (or VCs) are called escaping [E], 
but among them we count escapes in hyperbolic [H] orbits. For 
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Table 3. The past motion of dynamically new comets. The table is organized in the same manner as Table [2] In the case of mixed swarm the mean values or 
deciles of Q and q are given for the returning part of the VCs swarm, where the escape limit of 120 000 au was generally adopted with the exception of four 
objects mark with an asterisk (C/2001 CI and C/2004 X3) or two asterisks (C/2001 K5 and C/2003 Gl), where the escape limit of 140 000 au and 200 000 au 
was applied, respectively. 



# Comet 


9osc 


£osc 


Number of VCs 


£2prev 


^/prev 








%of 








au 


deg 


[Rl 


[El 


[HI 


10 3 au 


au 


lO^au^ 1 


dyn. new 
























10 au- 


15 au 


- 20 au 


[1] [2] 


[3] 


[4] 


[5] 


[6] 


[71 


[81 


[9] 


[10] 






[11] 




4 


C/1974 VI 


6.02 


-28.0 


2739° 


2262 


337 


54.1 - 80.6 - 112.7[R] 


12.49 - 78.23 -576.36[R] 


17.7±12.1 


96.9 


- 92.0 ■ 


88.8 


7 


C/1978 Al 


5.61 


31.4 


4814" 


187 


5 


37.8 - 52.7 - 82.2[R] 


9.01 - 19.52 - 141.32[R] 


37.5±11.9 


85.1 


- 63.3 ■ 


50.9 


8 


C/1978 G2 


6.28 


34.2 


792 


4209 


3585" 


33.3 - 62.5 - 109.0[R] 


7.25 - 27.05 - 568.1[R] 


-22.4±37.8 


95.7 


- 94.0 ■ 


93.1 


12 


C/1984 W2 NG 


4.00 


-33.0 


3488" 


1513 


43 


60.8 - 84.4 - 113.1[R] 


12.8 - 91.3 - 606.0[R] 


20.3± 


8.6 


95.1 


- 91.6 - 


88.0 


15 


C/1987 W3 


3.33 


-64.7 


4534" 


467 


1 


58.5 - 78.0 - 110.7[R] 


14.1 -61.1 -451.8[R] 


24.7± 


7.3 


95.8 


- 89.8 - 


82.7 


16 


C/1988 Bl 


5.03 


47.7 


3770" 


1231 


7 


66.7 - 89.6 - 116.9[R] 


39.9 - 208.3 - 985. 0[R] 


20.2± 


7.0 


99.9 


- 99.1 - 


98.1 


17 


C/1992 Jl 


3.00 


43.7 


5001 








73.6 ± 2.5 


20.4 - 29.7 - 43.9 


27.2± 


0.9 


100- 


99.0 - 


91.1 


20 


C/1997 Al 


3.16 


19.9 


4993 


8 





83.4 - 91.6 - 101. 2[R] 


56.7 - 111.2 - 224.6 [R] 


21. 8± 


1.7 


100 


- 100- 


100 


21 


C/1997 BA^ G 


3.44 


-27.7 


5001 








60.1 - 62.8 - 65.8 


15.9- 19.5-24.7 


31. 8± 


1.2 


100- 


95.3 - 


44.6 


26 


C/1999 J2 


7.11 


49.8 


5001 








90.4 ± 2.6 


160- 199-253 


22.1± 


0.6 


100 


- 100- 


100 


27 


C/1999 K5 


3.26 


-34.0 


5001 








93.6 ±4.1 


158 - 219 - 314 


21.4± 


0.9 


100 


- 100- 


100 


31 


C/1999 U4 


4.92 


30.7 


5001 








62.9 ± 1.0 


37.3 ± 2.7 


31. 8± 


0.5 


100 


- 100- 


100 


33 


C/2000 Al 


9.74 


33.2 


5001 








49.3 ± 2.4 


21.4-24.9-29.7 


40.6± 


2.0 


100 


-100- 


97.8 


39 


C/2001 CI* 


5.10 


8.3 


4019" 


982 





106.2 - 122.2 - 136. 9[R] 


68.9 - 179.0 - 366.0[R] 


15.9± 


2.1 


100 


- 100- 


100 


41 


C/2001 K3 


3.06 


-17.6 


4919" 


82 





50.2 - 64.0 - 87.0[R] 


6.08 - 14.9 - 83.7[R] 


31. 1± 


6.8 


67.8 


- 50.6 - 


40.2 


42 


C/2001 K5** 


5.18 


30.2 


4741" 


260 





182.0 ± 4.8 


19200 ± 2300 


9.6± 


0.4 


100 


- 100- 


100 


43 


C/2002 A3 


5.15 


9.0 


4942" 


59 





86.3 - 95.9 - 107.6 


48.5 - 82.9- 156.2 


20.7± 


1.8 


100 


- 100- 


100 


44 


C/2002 J4 


3.63 


-32.5 


5001 








58.8 ±2.4 


22.5 - 27.7 - 35.3 


34.1 ± 


1.4 


100 


-100- 


97.7 


46 


C/2002 L9 


7.03 


-51.7 


5001 








54.7 ± 1.3 


19.4-21.4-24.1 


36.5± 


0.9 


100 


-100- 


80.9 


48 


C/2003 Gl** 


4.92 


21.9 


5001 








136.6- 144.4- 152.9 


2041 - 2870 - 4056 


13.7± 


0.6 


100 


- 100- 


100 


49 


C/2003 S3 


8.13 


-11.5 


5001 








50.3 -55.5 -62.1 


14.0- 16.8-22.1 


36.0± 


2.9 


100- 


■ 78.0 - 


19.6 


51 


C/2004 PI 


6.01 


-20.2 


5001 








57.0-64.1 -73.3 


18.1 -27.8-51.0 


31. 1± 


3.0 


100- 


97.9 - 


82.4 


52 


C/2004 T3 


8.86 


-34.7 


5001 








39.1 -43.8-49.8 


13.4- 15.9-20.7 


45.7± 


4.2 


100- 


64.8 - 


12.5 


53 


C/2004 X3* 


4.40 


36.6 


2772° 


2229 





118.6- 134.0- 144.9[R] 


746 - 1677 - 2755[R] 


13.4± 


2.1 


100 


- 100- 


100 


54 


C/2005 B1 NG 


3.20 


17.8 





5001 









4.0± 


0.6 


100 


- 100- 


100 


56 


C/2005 Gl 


4.96 


55.3 


4551" 


450 





119.0±6.0[R] 


490 - 739 - 1046[R] 


16.6± 


1.0 


100 


- 100- 


100 


57 


C/2005 K1 NG 


3.69 


17.1 


112 


4889" 





106.5- 114.9- 118.8[R] 


207.2- 365.1 -465.2[R] 


11. 0± 


2.7 


100 


- 100- 


100 


58 


C/2005 Ql 


6.41 


8.9 


5001 








78.8-87.8-99.7 


34.6-56.0- 105.3 


22.7± 


2.0 


100 


- 100- 


100 


59 


C/2006E1 


6.04 


-31.5 


5001 








56.6-61.5 -67.6 


21.1 -29.3 -45.0 


32.5± 


2.2 


100- 


99.8 - 


93.3 


60 


C/2006 Kl 


4.43 


-64.0 


4621" 


376 





114.0- 121.7- 129.2[R] 


375 - 558 - 793[R] 


16.2± 


0.9 


100 


- 100- 


100 


64 


C/2007 Yl 


3.34 


35.6 


4631" 


370 


10 


39.2 - 56.8 - 90.0[R] 


4.6 - 14.5 - 230.0[R] 


34.1±1 


2.5 


63.3 


- 52.9 - 


47.3 



Table 4. The past motion of 7 comets with uncertain new/old status. The table is organized in the same manner as Table [2] In the case of mixed swarm the 
mean values or deciles of Q and q are given for the returning part of the VCs swarm. The escape limit of 1 20 000 au was used for all these comets. 



# Comet 


Rose- 




Number of VCs 


(2prev 


^/prev 


1/flori 


% of 




au 


deg 


[Rl 


[El 


[HI 


10 3 au 


au 


10- 6 au-' 


dyn. new 




















10 au - 15 au - 20 au 


[1] [2] 


[3] 


[4] 


[5] 


[6] 


[71 


[8] 


[9] 


[10] 


[11] 


19 C/1993 K1 


4.85 


-2.3 


4753" 


248 


2 


51.2- 68.3 -94.6[R] 6.33 - 10.14 - 32.89[R] 29.1± 7.7 


53.4-33.0-23.8 


23 C/1999 Fl 


5.79 


-15.2 


5001 








53.5 ± 1.0 


12.12 ±0.48 


37.4± 0.7 


100-0-0 


35 C/2000 Kl 


6.28 


21.5 


5001 








50.2 ± 2.9 


12.6- 15.0- 18.5 


40.0± 2.3 


100 - 49.3 - 4.8 


38 C/2000 Yl 


7.97 


-28.0 


5001 








30.4- 33.1 - 36.4 


9.8- 10.5- 11.6 


60.4± 4.2 


80.9 - - 


40 C/2001 Gl 


8.24 


46.8 


5001 








49.5 ± 3.4 


12.4- 14.5- 18.0 


40.6± 2.8 


100 - 40.7 - 4.4 


50 C/2003 WT 42 5.19 


57.0 


5001 








44.03 ±0.31 


12.45 ± 0.22 


45. 4± 0.3 


100-0-0 


62 C/2006 YC 


4.95 


27.7 


4965" 


36 





32.8 - 43.9 - 64.6[R] 


6.5 - 9.9 - 34.6[R] 


45.7±12.0 


49.9 - 28.7 - 20.0 



two comets in our sample, namely C/2001 CI and C/2004 X3, 
we decided to increase the threshold value up to 140 000 au be- 
cause the median of aphelia of their orbits were slightly below this 
value (both are marked with (*) in Table [3}. For the past motion of 
next two comets, C/2001 K5 and C/2003 Gl, we present the res- 
ults for the escape limit of 200 000 au, see Table [3] where they are 
marked with (**). For such a huge threshold value all VCs of these 



comets are returning and we are able to present more reliable de- 
scription of the previous perihelion distance then for the standard 
threshold, where we are restricted to the synchronous variant (see 
below) because all VCs are escaping. The classification of all these 
four comets as dynamically new is by no means influenced by these 
escape border extensions due to very large previous perihelion dis- 
tances of all of them. 
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Figure 6. Past and future heliocentric distance changes of all VCs repres- 
enting C/2005 Kl. 

For the detail description of the dynamical model as well as its 
numerical treatment the reader is kindly directed to Paper I. Basing 
on the conclusions from that work we used both, Galactic disk and 
Galactic centre terms in all calculations. Comparing with Paper I 
all the parameters of the Galactic gravity field are kept unchanged, 
including the local disk mass density, p = 0. 100 M,: /pc 3 .The rules 
of stopping the numerical integration were as follows: if all VCs for 
a particular comet were returning, all of them were stopped at indi- 
vidual previous/next perihelia. There were also a synchronous vari- 
ant, in which all VCs were stopped simultaneously when the nom- 
inal VC reached previous/next perihelion. When all VCs were es- 
caping the calculation was always terminated synchronously, when 
the fastest VC crossed the escape limit, usually equal to 120 000 au. 
If for a particular comet the swarm of VCs consists of both re- 
turning and escaping VCs, the returning part was stopped at previ- 
ous/next VC perihelia and the rest (escaping ones) when the fastest 
escaping VC crossed the escape limit. For these mixed swarms we 
also performed a synchronous variant, in which all VCs (both re- 
turning and escaping) were halted when the fastest VC crossed the 
escape limit. 

In Tables l2l3l and l4l in columns [8] and [9] we presented pre- 
vious aphelion and perihelion distances, respectively. In Table [9] 
the analogous data for the comets returning in the future are given 
in columns [7] and [8]. Depending on the distribution character- 
istics we use two different ways of aphelion and perihelion dis- 
tance presentation: if the distribution can be reliably approximated 
with the Gaussian one, we present the estimated mean value and its 
standard deviation. The example of Gaussian distribution of past 
elements one can find in Fig. [3] In the case of highly deformed dis- 
tribution, we present three deciles: 10th, the median and 90th. The 
example of non-Gaussian distribution of past elements one can find 
in Fig.H 

3.1 Overall statistics 

For the past motion we obtained 42 comets with all VCs return- 
ing, 19 comets with the mixed VC swarms and only 3 comets 
fully escaping (but with all VCs on highly eccentric elliptical or- 
bits). For statistics presented in this section the standard escape 
limit of 120 000 au for all investigated comets were adopted. Al- 
most all 19 comets with mixed VC swarms have the majority of 
returning clones, only two swarms consist mainly of escaping VCs 



Table 5. Overall VC distributions for previous and next perihelion pas- 
sage (based on the standard escape limit of 120 000 au) for all investigated 
comets except C/1999 H3, see text. 



deciles 


10 per cent 


50 per cent 


90 per cent 


'/previous 


4.11 au 


11.83au 


129.60 au 


^previous 


29000au 


48400au 


92 400 au 


time to previous perihelion 


9.66 Myr 


3.74 Myr 


1.82 Myr 


'/next 


3.00 au 


4.85 au 


8.30 au 


(2next 


3 370au 
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(C/2005 Kl and C/1978 G2, the last comet is the only one with 
nominal hyperbolic original orbit). 

In total, for the past motion of studied comets we obtained 
the 275 042 returning VCs (87.3 per cent) in 315 063 of all start- 
ing from the Oort spike (outside Oort spike was NG swarm of 
C/1999 H3). Statistics of previous perihelion and aphelion distri- 
butions for all these returning VCs are shown in Table [5] 

It should be stressed however, that we call escaping all VCs 
moving further than 120 000 au from the Sun. This is motivated 
by the fact, that because of a large heliocentric distance and huge 
orbital period these VCs may have their orbits modified by (even 
weak) stellar perturbations in the past (or future). The only possib- 
ility is to state, that at the moment their dynamical history is im- 
possible to be revealed. But one should not interpret these comets 
as of interstellar origin. The vast majority of the VCs escaping in 
the past still move in elliptical, heliocentric orbits and we have no 
direct evidence that they were not the Solar System members. See 
Section EOl for a detailed analysis of some particular examples. 

It is widely known, that the situation is quite different as it 
concerns the future motion. In Paper I the great majority of the 
investigated 22 comets with q < 3.0 au (about 77 per cent) were 
ejected from the Solar System by planetary perturbations. In the 
present sample of 64 comets this percentage is significantly smal- 
ler: 33 comets (about 52 per cent) are ejected in the future, two 
comets from Table|9](see below) and all from TablellOl 

We obtained 31 comets with all VCs escaping in the future 
on hyperbolic orbits (with the exception of a small part of the VC 
swarm of C/1987 HI escaping on the extremely eccentric elliptical 
orbits). We also obtained 27 comets fully returning in the future 
and only 6 comet with the mixed swarms. All these mixed swarms 
mainly consist of escaping VCs, nominal VC of two comets have 
hyperbolic future orbit. Thus, it seems probable that these two 
comets (C/1978 G2 and C/2006 S2, see Table[9} also are escaping 
from the solar system. In case of two comets with the mixed future 
swarm (but without any hyperbolic VCs, namely C/1997 J2 and 
C/2002 J5) it is possible to obtain a fully returning future swarm 
by applying the escape threshold enlarged up to 200000 au. For 
C/2002 J5 this is illustrated in Fig. [5] where we present the helio- 
centric distance changes of all VCs representing C/2002 15, both 
one orbital period to the past and future. In contrast to Fig|6]the past 
VCs swarm is very tight here and the future swarm, while crossing 
the standard escape border of 120000 au, is all returning, having 
the future perihelion distance greater than 2 500 au. 

For the future motion, we have in total 135 661 (43.1 per 
cent) returning VCs and their statistics is presented in lower part 
of TableU] 
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Figure 7. Past and future orbital evolution of the nominal VC for 
C/1976 D2, an example of Oort spike comet detectable at three consecutive 
perihelion passages. 



-3x10 8 



-2xl0 6 



Time [years] 



Figure 8. Past and future orbital evolution of the nominal VC for 
C/1980 El, see text for a detailed description. This comet passed previous 
perihelion 2.6 Myr ago at 2 au and it will be ejected in the future in hy- 
perbolic orbit. The effect of strong planetary perturbation is easily visible 



3.2 Comets with extremely large semimajor axes in the past 

In our sample, we have six comets with extremely large ori- 
ginal semimajor axes (l/a or j < 15 x 10 -6 au -1 ). Three of them, 
C/2001 K5, C/2003 Gl and C/2005 Bl, have well-determined 
1 /a or i and completely escaping swarms of VCs for the standard 
escape limit of 120 000 au used in this paper. The first two have 
similar osculating perihelion distance of about 5 au and a very sim- 
ilar behaviour in the past. It can be seen from Table [3]that a swarm 
of C/2003 Gl is completely returning for the escape limit shifted to 
200 000 au with the previous perihelion of about a few thousand au 
from the Sun. The same is true in the case of comet C/2001 K5 for 
slightly larger escape border but in this case the previous perihelion 
distance is almost one order greater. This rather suggest, that both 
are Solar System members, but due to their large semimajor axes 
and long orbital period we cannot exclude that their past motion 
was disturbed by stellar perturbations and as a result their dynam- 
ical history was quite different. 

Third comet with completely escaping past swarm, 
C/2005 Bl, has formally semimajor axis of about 
250 ± 50 thousand au! This comet is also unique in the sense 
that NG effects determined for this comet caused elonga- 
tion rather than shortening of its NG original semimajor axis 
relative to GR solution (see Section 12.11 and Fig. |2j, how- 
ever the difference between NG and GR models is small: 
A(l/a) = (l/fl) ri,NG - (V«)ori,GR = (-3.2 ±0.9) x lO^au" 1 . 
Taking all this into account comet C/2005 B 1 seems to be unique 
and we cannot even rule out its interstellar origin. It seems worth 
to mention that this comet is the only one of six comets considered 
in this section that is returning in the future. 



The remaining three comets with extremely large past se- 
mimajor axes are C/1978 G2, C/2004 X3 and C/2005 Kl. All three 
have mixed past swarms of VCs but while C/2005 Kl and C/2004 
X3 have all VCs in elliptical orbits, the majority of VCs repres- 
enting C/1978 G2 are hyperbolic. However it must be noted, that 
C/1978 G2 has poorly determined orbit. In fact it is the worst de- 
termined orbit throughout our sample, what comes from an ex- 
tremely small number of observations - we have only 7 positions 
of that comet. As a result, its past swarm of VCs is very dispersed. 
Nevertheless, it is the only one comet in our sample with l/a or i 



formally negative. Having 72% of past VCs (including a nominal 
one) moving on hyperbolic orbits this comet seems to be a candid- 
ate for an interstellar one. 

The past and future dynamical evolution of C/2005 Kl is 
shown in Fig. [6] We plot here the heliocentric distance of all 5001 
VCs with respect to the time. The zero point in the time axis corres- 
ponds to the observed perihelion passage of this comet (2005 Nov. 
21). To obtain such a plot for this particular comet, we allowed all 
VCs to move as far as 500 000 au from the Sun, what takes more 
than 80 million years for some of them. Even such an extremely 
distant escape limit is not sufficient to obtain purely returning past 
swarm of VCs for this comet. While all orbits are elliptical, their 
high dispersion forced us to conclude, that the dynamical history of 
C/2005 Kl cannot be determined basing on available observations. 
In contrast, all its future VCs are ejected from the Solar System 
in hyperbolic orbits (nominal l/flf ut = (— 82.4±3.3) x 10~ 6 au _1 ) 
without any doubt. There is an additional interesting detail in its 
past evolution depicted in Fig. [6] We marked with black lines the 
future motion of the nominal VC (a-curve), the past motion of the 
nominal VC (b-curve) and the past motion of one additional VC 
(c-curve), which represents an interesting dynamical scenario. Its 
orbital period is equal to the period of long term perihelion dis- 
tance changes due to the Galactic tides. As a result its previous 
perihelion distance can be arbitrarily small. However, because it 
takes some 57 million years to reach previous perihelion for this 
VC, one should treat this scenario as rather questionable due to po- 
tential stellar perturbations suffered by it during such a long time 
interval. 

The past evolution of VC swarm of C/2004 X3 is much more 
tight than that of comets C/1978 G2 and C/2005 Kl. However, to 
obtain the majority of clones coming back it was necessary to shift 
the escape limit to 140 000 au (see Table[3) and only unrealistically 
large escape limit of 260 000 au gives a whole VC swarm returning. 

3.3 Small previous perihelion passage distances 

In the light of Jupiter- Saturn barrier concept one should expect, 
that most of the Oort spike comets should have the previous peri- 
helion distance well out of rich of planetary perturbations. This is 
not the case. In the sample of 22 small perihelion comets (q osc < 
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Table 6. Comets with previous perihelion distance smaller than the observed one. This means that they were observed after the minimum point in the Galactic 
evolution of the perihelion distance, see also Fi g 1 1 5 1 In column 5 the statistics of perihelion distance changes of all VCs during the last orbital revolution is 
presented by three deciles, and the percentage of negative Aq in column 6. Since changes in q are highly correlated with the Galactic argument of perihelion 
CO, its evolution is also shown in columns 7 and 8, as well as the latitude of the perihelion direction b n. For comparison the previous perihelion distance in the 
galactic tide disk model only is given in parentheses in the third column 
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3.0 au) investigated in Paper I we obtained previous perihelion dis- 
tances g pl -e V < 15 au for 15 comets (almost 70 per cent). Now, 
in the sample of 64 LPCs with q osc > 3 au we obtained signific- 
antly smaller fraction of such comets, but still almost 50 per cent 
of the sample have the previous perihelion distance smaller than 
15 au. Moreover, among them, 6 comets (C/1972 LI, C/1976 D2, 
C/1976 Ul, C/1979 M3, C/1980 El and C/1997 J2) have the 
previous perihelion distance smaller than the osculating one (see 
Table |6) ! When comparing with Paper I, the percentage of comets 
observed at the greater perihelion distance than the previous one 
is roughly the same, but the ^-changes for small perihelion comets 
are smaller, below 0.3 au. It should be stressed here, that observing 
LP comets during the increasing phase of their perihelion distance 
evolution is the direct evidence that their dynamical history fol- 
lowed one of the two possibilities: either they were strongly per- 
turbed by planets during their previous perihelion passage (what 
switched the phase of the perihelion distance evolution) or they 
have moved unperturbed through the Jupiter-Saturn barrier in the 
past. Of these six comets, C/1976 D2 suffered practically no plan- 
etary perturbations in the observed perihelion passage (see Fig. [7J- 
These comets can also be found in the very bottom part of Fig. 1151 
with some theoretical inteipretation of this distribution given in 
Section|43l 

An example of such a situation is depicted in Fig. [8] where 
past and future dynamical evolution of C/1980 El is illustrated by 
the evolution of orbital elements of its nominal VC. Like for any 
other VC we followed numerically its motion under the influence 
of Galactic perturbation. Since we present several similar plots, we 
describe it here in more detail. For the past motion we started from 
the original orbit while for the future motion from future orbit, both 
were obtained with the NG effects included in this specific case, see 
Section [2~T1 for additional information. The horizontal axis shows 
the moment of osculation, for which orbital elements are calcu- 
lated and plotted, the zero point corresponds to the observed perihe- 
lion passage. The left vertical axis is expressed in au and describes 
both heliocentric distance of a VC (r, thin vertical lines) and its 
perihelion distance (q, continuous line). The right vertical axis de- 
scribes angular elements (calculated in the Galactic frame) and is 
expressed in degrees. We plot here the synchronous evolution of 
the argument of perihelion (to, dashed line) and inclination (z', dot- 
ted line). The thick lines depict the real dynamical VC evolution 
while their continuation with the thin lines depicts its potential mo- 
tion in the absence of planetary perturbations in the previous/next 
perihelion. All grey (cyan) lines (right from the zero point) describe 



additionally an artificial variant of the future motion, in the absence 
of all planetary perturbations during the observed perihelion pas- 
sage too. The discontinuities of the thick lines at the zero point of 
the time axis are the result of a close encounter of C/1980 El with 
Jupiter ( A = 0.228 au, 9.46 December 1980). 

The original perihelion distance of C/1980 El q or \ =3.17 au 
while the previous one (almost 2.6 million years ago) ^ prev = 
2.16au. It should be noted, that at previous perihelion passage 
C/1980 El was perturbed by planets, possibly rather strongly due 
to its small previous perihelion distance. Since it is impossible to 
calculate this perturbation (due to large uncertainties in planetary 
positions 2.6 million years ago) one should treat the orbital ele- 
ment evolution left from the previous perihelion passage likely to 
be completely fictitious so we pointed them in thin lines. 

The largest perihelion distance increase one can found for 
C/1976 D2, where the observed value was 6.88 au but the previous 
one 5.50 au (see Table|6](. This comet is also a clear evidence, that 
moving through the Jupiter- Saturn barrier can be observed. Planet- 
ary perturbations in this case were very weak, slightly decreasing 
its inverse semimajor axis from 1 /a or j = (56.9 ±7.3) x 10~ 6 au _1 
to l/af ut = (54. 1 ±7.3) x 10~ 6 au _1 . This almost unperturbed mo- 
tion through perihelion is depicted in Fig(7] In fact, this comet is 
a 'double evidence'. First, we observed it at larger perihelion dis- 
tance than the previous one, with all consequences described above. 
Second, the observed perihelion passage demonstrated an unper- 
turbed motion through the Solar System and its next perihelion dis- 
tance is even larger, while still observable! This is discussed in a 
more detail in the next section. 



3.4 Small next perihelion passage distances 

While the Jupiter-Saturn barrier mechanism predicts that great ma- 
jority of all LP comets that approach the Sun closer than 10-15 au 
should be definitely removed from this population, we observe that 
over 40 per cent of our sample ( 26 comets ) will keep moving on 
the typical LP comet orbits with small (q < 10 au) next perihelion 
passage distances. Moreover, six of them remain to be members 
of the Oort spike (l/a fut < 10~ 4 au). These comets (C/1976 D2, 
C/1999 Fl, C/2000 Al, C/2002 L9, C/2004 T3 and C/2005 Gl) 
constitute an important, direct evidence, that about 10 per cent 
of the large perihelion distance Oort spike comets can move dir- 
ectly through the Jupiter- Saturn barrier and remain observable. It 
is worth to mention that for a hypothetical observer of such comets 
next perihelion passage they can successfully pretend to be a res- 
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Figure 9. Past and future orbital evolution of the nominal VC for 
C/1999 Fl, see text for a detailed description. Another example of Oort 
spike comet moving through the Jupiter-Saturn barrier and remaining ob- 
servable LP comet. 50 per cent of VCs representing this comet will have 
the next perihelion distance smaller than 0.012 au and 90 per cent smaller 
then 0.27 au. 



Table 7. Observed perihelion distribution in the sample of large perihelion 
Oort spike comets (for explanation of dynamically new, dynamically old 
and dynamically uncertain comets see Section[4j. 
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Table 8. Observed inclination distribution in the sample of large perihelion 
Oort spike comets (for explanation of dynamically new, dynamically old 
and dynamically uncertain comets see Section[4j. 
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ult of the lKaib & Ouirml J2009h scenario discussed in Section 14.41 
Three of them have their semimajor axes significantly shortened, 
what makes the perihelion distance evolution under the Galactic 
tides much slower. 

An impressive example of very small next perihelion distance 
is the case of C/1999 Fl, see Fig. [9] The previous perihelion dis- 
tance of this comet was 12. 1 ±0.48 au, the observed one 5.79 au but 
the next with high certainty will be smaller than 0.3 au! In contrary 
to the most probable scenario attributed to the Jupiter-Saturn barrier 
crosser, the semimajor axis of this comet was slightly increased due 
to the planetary perturbations (the same happened to C/1976 D2) 
and as a result the next perihelion passage will be closer to the Sun 
than in the absence of planets. Typical planetary perturbation here 
are very small. In Paper I we obtained only 9 returning comets for 
the future motion (from the sample of 22 comets) but all of them 
have perihelia inside the observable zone. 



3.5 Comets with next semimajor axis below 2 000 au 

According to our analysis ~ 12 per cent (8 objects) of the observed 
large perihelion Oort spike comets return in orbits similar to that of 
comet C/1996 B2 Hyakutake ((l/a) fut = 554 x 10~ 6 au" 1 ), how- 
ever that comet had the original orbit more tightly bound than future 
orbit. Six of these comets creates a noticeable local maximum in 
the 1 /«f ut distribution displayed in Fig. [12] This maximum consists 
of two dynamically new comets, three dynamically old comets and 
one with an uncertain past dynamical status (see the definitions in 
Section l4~TT ). Three of these objects have g osc < 3.5 au. The shortest 
future orbit is for C/2002 A3, which future semimajor axis equals 
~ 162 au. This comet will return in about 1600 years with the peri- 
helion distance of 5.15 au. 



4 NEW AND OLD LP COMETS 

The terms of new and old long-period comets are widely used in 
literature for many decades. Sometimes the authors add adjectives: 
dynamically or physically, to inform the reader, what criteria they 
use to distinguish between new and old LPCs, but in most cases, 
the intention is that dynamically new should appear as physically 
new and vice verse. Historically, the first criterion used in this field 
was simply the semimajor axis, a, value. The widely accepted state- 
ment was, that all comets with a > 10 000 au were dynamically new. 
This was used for example by Oort jl950h when introducing his 
concept of the di stant comet reservoir, n ow called the Oort Cloud. 
Just a year later, lOort & Schmidj dl95ll) published a paper which 
seems to be the source of the widely quoted and repeated opinion, 
that new LPCs are more active and brighter. In fact, nowadays it is 
very difficul t to pr ove the truth of this statement, see for example 
IDvbczvnskil feOOlh . who collected large number of LPCs absolute 
magnitudes and fo und no correlation with their dynamical history. 
IDvbczvnskil ( 1200 lh also showed, that using a > 10 000 au as the cri- 
terion of being the dynamically new LPC seems to be completely 

unsatisfactory. 

Recently, iFinkl fe009h presented an extended taxonomic sur- 
vey of comet composition, based on their spectroscopic observa- 
tions. As it concerns LP comets, he also did not found any correla- 
tion with the semimajor axis (see for example Fig. 8 in the quoted 
paper). 

4.1 How we distinguish them? 

Let's start with definitions. We use the term dynamically old 
LP comet for objects with previous perihelion passage distance 
smaller than some threshold. This threshold value should describe 
the sphere of significant planetary perturbations. In Paper I we used 
15 au as this limit but now we decided to use three different values, 
namely 10, 15 and 20 au in parallel to observe how the new/old 
classification depends on it for investigated comets. Because we 
replaced each individual comet with a swarm of 5001 VCs we ap- 
plied the above-mentioned criterion individually to each VC and 
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Figure 10. Past and future orbital evolution of the nominal VC for 
C/1999 J2. An example of dynamically new comet; from the swarm of its 
VCs it comes, that 90 per cent of them had their previous perihelion distance 
greater than 160 au. 
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Figure 11. This Aitoff projection sky map shows the distribution of 
aphelion directions for all (discovered before 2009) large perihelion dis- 
tance Oort spike comets in galactic coordinates. Squares, circles and tri- 
angles show comets investigated in this paper where symbols indicate dy- 
namically old comets, dynamically new and comets of uncertain dynamical 
status, respectively; full and open markers represent comets returning and 
escaping in the future. Black crosses show 11 comets discovered before 
1970. 



then classified a particular comet depending on the percentage of 
escaping VCs - if more than 50 per cent of VCs were escaping in 
the past we call the parent comet a dynamically new one with re- 
spect to the particular threshold value. It is worth to mention, that 
except of a few cases this percentage is significantly closer to zero 
or 100 per cent, see col. [11] in Tables[2]-S With this definition, 
a dynamically new LP comet should have moved (before the ob- 
served perihelion passage) in the orbit which is free from planetary 
perturbations and therefore can be used to study the source region 
of LP comets by tracing its motion back in time under the Galactic 
perturbations. From the point of view of the above-mentioned three 
different threshold values and basing on their motion in the past we 
finally divided the whole sample of 64 comets into three groups. 
In the first one, see Table [2] we placed 26 comets dynamically old 
with respect to all three values. The second group, see Table(3] con- 
sists of 3 1 dynamically new comets with respect to all three values 
(26 comets) or only to the two lower values (10 au and 15 au; 5 



comets.) The third one, see Table [4] groups 7 comets of the uncer- 
tain dynamical age: they are new if one take 10 au as the threshold 
value, but old for the greater threshold values. 

Observed perihelion distance and ecliptic inclination distribu- 
tions of all investigated comets one can find in the first rows of 
Tables [7}{8] There is a clear observational selection signature in 
the perihelion distance distribution. The ecliptic inclination distri- 
bution presented in Table [8] is shown for two groups of comets: 
with q < 4.5 au and q ^4.5 au. In the whole sample of the large 
perihelion Oort spike comets discovered since 1970 there are more 
comets moving in prograde orbits (56.2 per cent) than moving in 
retrograde orbits (43.8 per cent), however, for the sub-sample of 
comets with 3.0 au^C q <4.5 au we observe the same number of pro- 
grade and retrograde orbit comets. It means that this disproportion 
comes entirely from a subset of comets with q^4.5 au, where the 
ratio of comets in prograde orbits to comets in retrograde orbits is 
about 1.6. (see Table[8). 

We have also found an interesting feature in the distributions 
presented in Table[8] While the number of prograde and retrograde 
orbits of comets with 3.0 < q < 4.5 au is equal, the proportion 
of dynamically new to dynamically old ones is reversed in these 
groups. In other words on prograde orbits there are about twice as 
many dynamically new comets than on retrograde orbits (9 : 4) and 
the proportion is opposite for the dynamically old comets (5:10). 

A nice example of comet dynamically new for sure, C/1999 J2, 
is presented in Fig.[l0] The previous perihelion passage of nominal 
orbit of this comet happened some 9 million years ago at the helio- 
centric distance of 200 au (80 per cent of its VCs have a previous 
perihelion distance in the interval between 160 and ~250au, see 
Table [3). 

The spatial distribution of aphelion directions of all Oort spike 
comets with q > 3 au (discovered before 2009) is presented in 
Fig. [TTJ Circles mark dynamically new comets, squares dynam- 
ically old ones and triangles show seven comets with the uncer- 
tain dynamical status, according to the definitions adopted in this 
section. In addition to these we included here 1 1 large perihelion 
distance comets (marked with crosses) discovered before 1970 and 
omitted in the present investigation. Such a spatial distribution is 
often used when searching for a signa ture of some specific per- 
turbers, see iMatese & WMtrmrel d201ll) and iFernandea d201 ll) for 
recent discussions of this subject. When looking for effects of a 
massive perturber moving on a distant heliocentric orbit one should 
expect a concentration of the aphelia directions along some great 
circle in the sphere. It seems to be difficult to interpret the distribu- 
tion presented in Fig.QTJin such a way. 



4.2 How are they different from each other? 

Figure [T2l shows the distribution of original and future 1/a as 
well as the distribution of planetary perturbations acting on comets 
during their passage through the inner solar system (A(l/a) = 
1/flfut — l/ fl ori)- The black and white parts of histograms repres- 
ent the dynamically new and dynamically old comets, and grey 
- comets with uncertain dynamical status in the sense described 
above (Section 14. lb . Firstly, we focus on the total distributions 
of cometary energies (left-hand side panels). It is clear that all 
three distributions visible in Fig.ll2l(a) show some deviations from 
the Gaussian model. However, the Gaussian fitting to the sample 
of 62 comets (C/1978 G2 and C/1999 H3 were excluded) give 
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Ml/a) = (l/a)„, - (l/a)„„ [AU-] A(l/a) = (l/a)„, - (l/a)„„ [AU-] A(l/a) = (l/a) w - (l/a)„„ [AU-»] 

(a) fee > 3.0 au (b) 3.0 ^ g OS c < 4.5 au (c) g osc ^ 4.5 au 



Figure 12. Distributions of original and future cometary energies measured by 1 / a 0[l (top panels) and 1 /flf ut (middle panels) for the observed sample of large 
perihelion Oort spike comets. The bottom panels represent the distribution of planetary perturbations acting on the comets during their passage through the 
planetary system (A(l/a) = 1 /<3f ut — l/a or i). The filled parts of histograms represent the dynamically old comets, and grey - the seven comets with uncertain 
dynamical status, i.e. comets not dynamically new for the limit of 15 au for previous perihelion distance but being dynamically new for the limit of 10 au. 
Dotted histograms in the middle and right-hand side panels ((b) and (c)) represent the total histograms of 64 comets from the left-hand side panel (a). 



10 100 1000 



0.0001 



D 

< 



10" 



T I I 1 1 1 1 1 I I I I ll| 1 1 I I I I llj 1 1 I I I I ll| 1 1 I I I I llj 

~«C/1999 H3 




Aq 



iprevious T-ori 



<Vi [AU] 



Figure 13. Aq vs l/a ml in logarithmic scales, derived for the observed 
sample of large perihelion Oort spike comets. The open dots show the dy- 
namically new comets, the filled dots - the dynamically old comets and grey 
squares - the seven comets with uncertain dynamical status, i.e. comets not 
dynamically new for the limit of 15 au for previous perihelion distance but 
being dynamically new for the limit of 10 au. The error bars for perihelion 
distances are described by deciles within 10-90 per cent, the vertical er- 
ror bars are described by 1-sigma error for normal distribution of l/a ri- 
Straight line represents the best fit to all presented points, except the Comet 
C/2006 S2. 

< l/a ori >= (37.8 ± 17.7) x 10 _6 au _1 with negative kurtosis3 
equal to -0.55. Generally 1 /a orl -distribution has a wider and lower 



We use standard definition of kurtosis: K = 
central moment, o~ is the standard deviation. 



- 3 , where 114 is the fourth 



peak around the mean value than normally distributed variable. 
Outside the horizontal scales of the middle and lowest panels are 
two comets that have suffered large planetary perturbations during 
their passage through the inner solar system: C/1980 El Bowell 
(A(l/a) = —16064 x 10 au _I mainly due to Jupiter encounter 
within 0.228 au on December 1980) and C/2002 A3 LINEAR 
(A(l/a) = +6153 x 10~ 6 au~' mainly due to Jupiter encounter 
within 0.502 au on January 2003). The planetary perturbations act- 
ing on the sample of these large perihelion comets show clear 
asymmetry relative to zero (see lower panel) with the negative me- 
dian value of A(l/a) = — 51.8 x 10~ 6 au~' whereas the distribution 
of 1 /flfut are more symmetric relative to zero (the middle panel) 
with median value of —6 x 10~ 6 au~'. Statistical analysis shows 
that the A(l/a)-distribution is closest to the Gaussian distribution. 
We estimated the value of mean planetary perturbations (represen- 
ted by the standard deviation of the A(l/a)-distribution) equal to 
285 x 10~ 6 au~' by fitting to the Gaussian distribution, however 
A(l/fl)-distribution has the non-zero kurtosis (0.201) and is asym- 
metric with the longer right tail (skewness equal to 0.237). The ob- 
tained mean planetary perturbation in comet energy is si gnificantly 
smaller than predicted by the numerical simulations ( Fernandea 
1 198 it iDuncanet aT1ll987r> . This is probably the result of the non- 
uniform inclination distribution in the observed sample of large 
perihelion comets (see Table [8), in contrast to the quoted simu- 
lations. With relatively large number of high inclination and retro- 
grade orbits the mean energy change is expected to be smaller. 

Future 1 /a-distribution seems to have the second small max- 
imum in the interval of 500 x 10~ 6 au _1 < l/flf ut < 600 x 
10~ 6 au -1 (middle panel of Fig.ll2t. This maximum consists of six 
comets (C/1974 Fl, C/1974 VI, C/1992 Jl, C/1993 Kl, C/1999 Ul 
and C/2000 CT54) where three of them have o osc < 3.5 au. 

Separate 1 /a-distributions of comets with moderately large q 
(3.0 ^ q OS c < 4.5 au) and very large q (q osc > 4.5 au) are shown in 
panels (b) and (c) of Fi gure 1121 where the contributions of all three 
dynamical groups of comets to original and future 1/a-distributions 
as well as to A(l/a)-distributions are presented. Both distribu- 
tions of original 1 /a are noticeably different. l/fl or i -distribution 
of comets with very large q is more compact with a clear max- 
imum around 1 /a 0[l ~ 40 — 50 x 10 au~' whereas the distribu- 
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Table 9. The future distributions of the returning and mixed swarms of VCs in terms of returning [R], escaping [E], including hyperbolic [H] VC numbers. 
Aphelion and perihelion distances are described either by a mean value for the normal distributions, or three deciles at 10, 50 (i.e. median), and 90 per cent. 
In the case of mixed swarm the mean values or deciles of Q and q are given for the returning part of the VCs swarm, where the es cape limit of 1 20 000 au was 
generally adopted with one exception of comet C/2002 J5, marked with an asterisk, where the escape limit of 140 000 au was applied. The upper-a index in 
columns 4-5 means that this part of mixed swarm includes the nominal orbit. Last column presents the value of 1 /af ut . 
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tion of comets with moderately large q is more flattened and is loc- 
ated somewhere between 30 — 40 x 10 au -1 . An additional dif- 
ference is reflected in the absence of very large perihelion comets 
with 1/flori > 70 x 10 au~' when we have four such comets 
with moderate q. It may be noted that future l/a or j -distribution 
and A(l/a or i)-distribution for moderate q comets are much more 
flattened while for comets with very large q we observe a clear 
maxima more or less around zero. Thus, we observe the excess 
of comets with very large q experiencing small (| A(l/a or j) |< 
10~ 4 au~') planetary perturbations. 

Figure Q/3] shows the relation between the changes in perihe- 
lion distance during the last orbital revolution and original semima- 
jor axes. Filled dots represent dynamically old comets, open dots 
- dynamically new comets and grey squares - seven comets with 
uncertain dynamical status. Six dynamically new comets with neg- 
ative Aq = g prev — q or { and four dynamically new comets with es- 
caping or almost escaping swarm are not included in this figure (see 
also Section [372t except comet C/2005 Kl which is shown as lowest 
point in the figure. Point standing off to the left of the fit line repres- 
ents comet C/2006 S2. It is worth noting that the comet C/1978 G2 



with a formally negative 1 /a or [ but large uncertainty of 1 /a „ -value 
is consistent with the presented fit (the logarithmic scale in the fig- 
ure makes it impossible to show that). Straight line represents the 
best fit to 53 points where the relation Aq ~ (l/a o ri)~ 4 ' 06± °' 16 was 
derived. Th e obtained expon ent is significantly smaller than that 
presented in lYabushital ( fl989h but remarkably clos er to the ex pect- 
ations of the first order Galactic disk tide theory, d&yl|[l_98^). It is 
still a little bit greater than theoretical but this might be the result 
of including the Galactic centre tide in our model. 



4.3 Evolution of orbital elements in Galactic frame 

It is a well known fact that under the separate Galactic disk tide the 
secular evolution of cometary orbital elements is strictly periodic 
and synchronous, i.e. the minimum of the perihelion distance coin- 
cides with the minimum of the Galactic inclination, with the max- 
imum of the eccentricity and with Galactic argument of perihelion 
crossing 90 or 270 degrees. Including of the Galactic centre term 
introduces only a small discrepancy from this regular patterns, but 
even during one orbital period in can manifest itself in some spe- 
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Table 10. The future distributions of comets with fully hyperbolic [H] swarms of VCs. Comets with NG effects are indicated by upper-NG index located behind 
the comet designation (column 2). In column 6 the eccentricity distribution at 1 20 000 au are described either by a mean value for the normal distributions, or 
three deciles at 10, 50 and 90 per cent. Last column presents the value of l/tff ut . 
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Figure 14. Galactic vs galactic latitude for previous perihelion passage or previous escaping on 120 000 au from the Sun (left-side and right-side pictures 
for dynamically new and dynamically old comets, respectively) or for the next perihelion passage or next escaping on 120 000 au from the Sun (middle picture 
for dynamically new comets). Right picture: almost all comets have galactic (O inside the first and third quarter; in the second and fourth quarter are five of six 
comets with q on > q ptm : C/1972 LI, C/1976 D2, C/1980 El, C/1997 J2, and C/1979 M3. 
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Figure 15. Galactic argument of perihelion co vs 1 / a or j (left panel) and vs Aq = qp rev — q ar [ (right panel) for all investigated comets (note the change from 
the linear scale to a logarithmic one in the middle of the vertical axis of the right panel). Circles and squares represent dynamically new and dynamically old 
comets, respectively, triangles - comets with uncertain dynamical status. Filled and open symbols mean returning or escaping in the future. Six comets with 
the negative Aq, i.e. with q m [ > qp av (see also Table[6}, can be located in the bottom part of this figure. These are (from the left to right side): C/1976 Ul, 
C/1972 LI, C/1976 D2, C/1980 El, C/1997 J2, and C/1979 M3. 



cific cases. This is true especially for small (Galactic) inclination 
orbits and for large orbital periods. 

Analysing the Galactic evo lution of cometary perihelion dis- 
tances [Matese&XIii^ii !15o4) introduced so called tidal charac- 
teristic '5' which tell us whether the perihelion distance of a partic- 
ular comet decreases (5 = — 1 ) or increases (5 = + 1 ) during the ob- 
served perihelion passage. In the dynamical model restricted to the 
Galactic disk tide 5 = — 5!g/i(sin(2co)). Due to the relatively short 
time intervals, limited to one orbital period we, despite of includ- 
ing the Galactic centre term in our calculations, do not observe any 
departure from this simple relation. It means, that for all comets in- 
vestigated here their perihelion distance decreased due to the com- 
bined Galactic tides when Galactic argument of perihelion CO is 
in the first or third quarter a nd increases oth erwise. As discussed 



Matese & Lissauerl d2004l) and recently in iMatese & Whitmirel 
J20T l|) this is the reason that dynamically new comets should have 
5 = — 1 (and ft) is in the first or third quarter) much more frequently. 
This, together with the distribution of the Galactic latitude of peri- 
helion direction b is illustrated in Fig.l 141 

As in the previous plots, circles denotes here dynamically new 
comets and squares dynamically old ones; in the left and middle 
panel filled symbols mark comets returning in the future while open 
symbols mark comets ejected in the future from the solar system 
as moving along a hyperbolic orbits. Due to the relation sinfo = 
sin ft) sin i comets may appear in this plots only between horizontal 
fe-zero axis and the solid, "triangle" lines. All angles are measured 
in the Galactic frame. 

In the left panel of Fig. [14] we present perihelion latitude 
verse argument of perihelion for nominal orbits of dynamically new 
comets at their previous perihelion passage. In accordance with the 
patterns described above all circles are located along solid borders 
of the allowed region, and only in the first and third quarter of CO 



what corresponds to the perihelion distance decreasing phase of its 
Galactic evolution (5 = —1). 

The same comets but at the next perihelion passage (i.e. two 
orbital revolution later) are plotted in the middle panel of Fig. Q5] 
One can observe, that almost all comets moved significantly and all 
four quarters are populated more or less uniformly. There is also an 
intriguing asymmetry in the displacements of escaping and return- 
ing comets. Please note, that the displacement of comets between 
the left and the middle panel incotporates both Galactic orbital 
evolution during two consecutive revolutions and planetary per- 
turbations during the observed perihelion passage. 

The right panel of Fig.[l4]is to be compared with the middle 
one. We present here angular elements of all (both returning and 
escaping) dynamically old comets at their previous perihelion pas- 
sage. If our definition of being dynamically old or new is relevant, 
new comets at their next perihelion ( only the filled points in the 
middle panel ) should have angular elements distribution qualit- 
atively similar to that of the old comets at their original perihelion 
passage (right panel). As one can easily observe we obtained a high 
level of such a similarity. 

IMatese & Lissauerl d2004h suggested that it should be a signi- 
ficant correlation between the tidal characteristics S ( i.e. the quarter 
of ft)) and the reciprocal of the semimajor axis 1 /a 0I i. As presented 
in Fig. [15] much more clear correlation can be observed between 
co and the change in the perihelion distance during the last orbital 
revolution: Aq = q prev — q OI [. It can be noticed that there are two 
distinct group of comets that can be found in all four quarters of 
co: dynamically new comets (that have largest Aq) and dynamically 
old comets with negative Aq (listed with details in Table |6j. The 
rest of the dynamically old comets as well as all dynamically un- 
certain ones can be found only in the first and the third quarter of 
co. 
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Table 11. The comparison of the previous perihelion distance (all values 
are expressed in au) of C/1984 W2, C/1993 Kl, C/1997 Al and C/1997 J2 
obtained in the present paper and in Dvbczvriski 1 2006). Demonstrated is 
both the influence of stellar perturbations and the Galactic centre tide term. 
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Figure 16. Past and future orbital evolution of the nominal VC for 
C/2001 K3. This comet came to the observability zone having a or j ~ 
32 000 au and might be considered as an example of the pr obable output 
of the mechanism proposed recently bv lKaib & Ouinnl ( 2009). 



4.4 New interpretations 

Recently an interesting ne w scenario of the LPCs dynamical evol- 
ution was pointed out by iKaib & Ouinnl fcOOSh . They showed, 
that in addition to preventing some LPCs to be observed, the 
Jupiter-Saturn barrier can help some inner Oort cloud comets (a < 
10000 au) to reach the observability zone. During the numerical 
simulation of the Oort cloud dynamical evolution they observed, 
that large percentage of the inner cloud comets follows a common 
pattern: due to weak Galactic perturbations their perihelia slowly 
drift towards the Sun and at heliocentric distances of 15-20 au they 
received a series of energy kicks from the giant planets increasing 
significantly their semimajor axis, see Fig. 1 in the quoted paper. At 
the late stage of this process a comet with q ~ 12 au goes outside 
the Planetary System along the orbit with a ~ 28000au and then 
it returns to the Sun at any small perihelion distance due to strong 
Galactic perturbations received during this last orbital revolution. 

They concluded that the majority of the observed LP comets 
can be produced by such a mechanism bu t except mentioning two 
objects in s trange orbits, 90 377 Sedna llBrown et al 1 120041) and 
2006 SQ372 dKaib et all20O9h . they do not provide any real comet- 
ary examples. 

In our study we obtained several results which can fit their 
prediction perfectly. By quick examination of Table |5]one can find 
that C/1978 Al, C/1997 BA 6 , C/2001 K3, C/2003 S3, C/2004 T3, 
C/2007 Yl are good candidates, as well as all 7 comets from 
Table |4] Previous perihelia of all these comets (as well as many 
VCs of other comets) lie in the vicinity of the Jupiter-Saturn barrier 
and it seems to be quite achievable that s ome of these comets w ere 
produced by the mechanism proposed bv lKaib & Ouin 

In Fig. [16] we present the past and future Galactic evolution 
of the nominal VCs of C/2001 K3. Its previous perihelion distance 
and original semim ajor ax is fits perfectly to the scenario proposed 
by Kaib and Quinn Visiting the planetary system 5.7 mil- 

lion years ago at the heliocentric distance of 1 1 au it probably re- 
ceived significant perturbations from Jupiter and Saturn. Therefore 
the probability that it was just placed on the large semimajor orbit 
being previously the inner cloud member seems to be quite signific- 
ant. In the future this comet was captured into small semimajor axis 
orbit of ~ 1 500 au, leaving permanently the Oort spike and having 
the next perihelion distance of equal to the observed one (3.06 au). 
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5 DISCUSSION AND CONCLUSIONS 

As a continuation of Paper I we attempted to characterize past 
and future motion of the next sample of the Oort spike (l/a or i< 
1 x 10~ 4 au -1 ) comets. In order to minimize possible biases due to 
indeterminable NG effects we decided to study here only comets 
having perihelion distance q osc > 3 au and precisely determined 
orbits. To this aim we omitted both 1 1 LP comets with g osc > 3 au 
but discovered before 1970 and three comets still (at the moment 
of this writing) potentially observable. Such a complete sample of 
64 large perihelion distance Oort spike comets allows us to obtain 
some statistical characteristics of this sample as well as individual 
past and future dynamics for all of them. In the process of osculat- 
ing orbit determination we succeed with NG effects detection for 
15 comets. Having a homogeneously obtained set of osculating or- 
bits we followed numerically their motion among planets back and 
forth, up to the heliocentric distance of 250 au, obtaining original 
and future orbits. Instead of integrating one orbit per comet we re- 
placed each body with a set of 5001 VCs and followed their motion 
individually. This allowed us to obtain all parameters of the original 
and future orbits together with their uncertainties. 



Then we analysed past and future motion of all VCs for one or- 
bital period, including both Galactic disk and Galactic centre tides 
and omittin g the perturbations from all known stars. The later was 
justified bv lDvbczynskil ( l2006h . who analysed by means of the ex- 
act numerical integration the influence of all known stars on the 
population of LPCs. In Table 6 of the quoted paper he listed 22 
long period comets for which stellar perturbation changed the pre- 
vious perihelion distance by more than 10 per cent. Only four of the 
large perihelion distance LPCs studied in the present paper can be 
found in that table. Thus in Tablell llwe can compare those results 
with the current calculations. There are several i mportant differ- 
ences between calculations of iDvbczvnskil j2006l) and the present 
paper. First, we included here th e Galactic centre ti dal term, not- 
ing its importance in some cases jDvbczvnskil d2006t) used only the 
disc tidal term). Second, we homogeneously determined all comet- 
ary orbits directly from observations and additionally included NG 
effects where possible. Third, we replaced each comet with the 
swarm of 5001 VCs, all compatible with the observations. This al- 
lowed us to observe the influence of the propagated observational 
uncertainties on the final results. 

In Table [TT] we present the current and previous results for 
four comets: C/1984 W2, C/1993 Kl, C/1997 Al and C/1997 J2, 
two of them have detectable NG effects (marked by the NG super- 
script after the comet designation). In the second column we in- 
cluded the perihelion distance value for the original orbits obtained 
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in this paper. Next two columns show the previous perihelion dis- 
tance for the nominal orbit of these comets obtained without (qj) 
and with (qdc) the Galactic centre tidal term. The next column of 
Table 1111. q prev , presents the result of the observational uncertain- 
ties propagated back to the previous perihelion. Last t wo columns 
show the previous perihelion distances obtained by IDvbczvhskil 
d2006l) : q\ presents the value obtained only using galactic disk tide 
while ^denotes the value obtained when perturbations of 21 most 
important stellar perturbers were included. One can easily note, that 
the VCs previous perihelion dispersion (column 5) is of one order 
greater that the differences between current and previous results, 
except in the case of comet C/1997 J2, where its current and past 
orbit can be determined with the great accuracy. It can be observed 
that changes in the dynamical model of the Galactic tides result 
in comparable (typically greater) changes in the previous perihe- 
lion distance than after the incorporating of the stellar perturbations 
(despite the fact, that these comets are the most sensitive for known 
stellar perturbations, so these differences should be treated as the 
extreme disparities). 

It should be noted, that even after including the action of all 
21 strongest stellar perturbers (still far to week to manifest them- 
selves) our classification of these four comets do not change in any 
way: C/1984 W2 and C/1997 Al are evidently dynamically new 
comets, C/1997 J2 is evidently dynamically old one and C/1993 Kl 
remains the uncertain case. As it concerns the comple teness of the 
stellar perturbers search the reader is kindly directed to IDvbczvhskil 
J2006h . where possible sources of this incompleteness were dis- 
cussed. We only summarise here these arguments stating, that the 
omission of any important (i.e. massive and/or slow and/or trav- 
elling very close to the Sun) star seems rather improbable but not 
impossible. If such a star will be discovered its dynamical influence 
on the long term dynamical evolution of all LPCs should be care- 
fully studied. This would also may change some particular results 
and conclusions presented here if it happened. 

Additionally it seems worth to mention that the comet 
C/1997 J2 belongs to an interesting group of six comets listed in 
Table [6] for which we obtained smaller previous perihelion dis- 
tance than the observed one. As it is clear from Table ITTI none of 
the compared dynamical models do not change this classification 
also. 

The swarms of VC stopped at 250 au from the Sun were next 
followed numerically for one orbital revolution to the past and fu- 
ture and as a result we obtained orbital characteristics of these ob- 
jects at the previous and next perihelion passages with respect to the 
observed one. Both Galactic disk and Galactic centre tides were 
taken into account for all comets. Basing on the previous perihe- 
lion distance we divided a sample of 64 large perihelion Oort spike 
comets into three groups: 26 dynamically old comets, 31 dynam- 
ically new comets and 7 comets with the uncertain dynamical age. 
When analysing orbits at the next perihelion we found that only 3 1 
comets will remain Solar System members in the future. The rest 
of our sample, 33 comets, will leave our planetary system along 
hyperbolic orbits due to the planetary perturbations. 

Detailed results and plots for all individual comets studied in 
this paper as well as summarizing cata logues of orbits will success- 
ively a ppear at this project web page, iDvbczvnski & Krolikowska 
feOllh . 

As a result of studying all individual dynamical evolutions as 
well as observing several interesting statistical characteristics of 
this sample of comets we can draw following main conclusions: 



(i) Observation selection and weighting are crucial for the pre- 
cise orbit determination. 

(ii) In contrary to popular opinion it is possible to determine 
orbits with NG parameters for some 25 per cent of large perihelion 
distance Oort spike comets. 

(iii) Incorporating these NG orbits makes the overall energy dis- 
tribution of our sample significantly different, modifying the shape 
of the Oort spike and position of the 1 /a „-peak. 

(iv) Replacing each individual comet with a swarm of virtual 
comets "equally well" representing the observations is a powerful 
method for analysing all uncertainties in the past and future mo- 
tion of LP comets. This allowed us for example to take the energy 
uncertainties into account when preparing the energy distribution. 

(v) In the absence of the recognized recent stellar perturbations 
it is quite possible to calculate previous perihelion distances of the 
LP comets, taking into account full Galactic potential. 

(vi) We have obtained clear correlation between the calculated 
change in the perihelion distance and the reciprocal of original se- 
mimajor axis. In contrast to some previous estimations we obtained 
the exponent quite similar to the theoretically predicted. 

(vii) On the basis of the obtained previous perihelion distance 
distributions we are able to distinguish comets that were observed 
for the first time, i.e. dynamically new, from the rest of comets vis- 
iting the observational zone at least twice - these are called dynam- 
ically old. 

(viii) By defining three different threshold values for strong 
planetary perturbations we can easily divide all 64 comets into 26 
dynamically new comets, 31 dynamically old comets and remain- 
ing 7 comets for which their dynamical status seems to be uncer- 
tain. 

(ix) With the analysis of the Galactic angular orbital elements 
evolution we found several significant fingerprints of the Galactic 
tides as a dominating agent delivering observable Oort spike 
comets at present. 

(x) In contrast to the overall picture of the so called Jupiter- 
Satum barrier we found that almost 50 per cent of our sample have 
the previous perihelion distance below 15 au. Moreover, we found 
several examples of comets which have moved through the Jupiter- 
Saturn barrier almost unperturbed. For six comets we found, that 
the observed perihelion distance was even larger than the previous 
one. 

(xi) Among future orbits we found 27 comets that will be ob- 
servable during their next perihelion passage. In contrast, 33 comets 
will be lost due to the hyperbolic ejection. 

(xii) We have also found, that among 64 large perihelion dis- 
tance Oort spike comets almost 25 per cent can be treated as pos- 
sible results of the recently proposed bv lKaib & Ouinnl J2009I) new 
source pathway from the inner Oort cloud to the observable zone. 



Acknowledgments. The research described here was partially sup- 
ported by Polish Ministry of Science and Higher Education funds, 
(years 2008-2011, grant no. N N203 392734). Part of the calcula- 
tion was performed using the numerical orbital package developed 
by Professor Grzegorz Sitarski and the Solar System Dynamics and 
Planetology group at SRC. The authors wish also to thank Pro- 
fessor Slawomir Breiter for valuable discussions on some partic- 
ular dynamical aspects of this research. We also thank Professor 
Giovanni Valsecchi who, as a referee, provided valuable comments 
and suggestions. This manuscript was partially prepared with LyX, 
the open source front-end to the TjX system. 



References 



Brown M. E., Trujillo C, Rabinowitz D., 2004, ApJ, 617, 645 
Byl J., 1986, Earth Moon and Planets, 36, 263 
Delsemme A. H„ 1987, A&A, 187, 913 

Dones L., Weissman P. R., Levison H. E, Duncan M. J., 2004, in: 

Comets II: Oort cloud formation and dynamics, pp 153-174 
Duncan M. J., 2009, Science, 325, 1211 
Duncan M. J., Quinn T., Tremaine S., 1987, AI, 94, 1330 
Dybczynski P. A., 2001, A&A, 375, 643 
Dybczynski P. A., 2002, A&A, 383, 1049 
Dybczynski P. A., 2006, A&A, 449, 1233 

Dybczynski P. A., Krolikowska M., 2011, 

http://apollo.as tro.amu.edu.pl/WCP 
Emel'yanenko V. V., Asher D. J., Bailey M. E., 2007, MNRAS, 

381, 779 

Fernandez J. A., 1981, A&A, 96, 26 

Fernandez J. A., 1994, in A. Milani, M. di Martino, & A. Cellino 
ed., Asteroids, Comets, Meteors 1993 Vol. 160 of IAU Sym- 
posium, Dynamics of Comets: Recent Developments and New 
Challenges, pp 223-+ 

Fernandez J. A., ed. 2005, Comets - Nature, Dynamics, Origin 
and their Cosmological Relevance Vol. 328 of Astrophysics and 
Space Science Library 

Fernandez J. A., 2011, ApJ, 726, 33 

Festou M. C, Rickman H., West R. M., 1993, Astronomy and 

Astrophysics Review, 5, 37 
FinkU., 2009, Icarus, 201, 311 

Kaib N. A., Becker A. C, Jones R. L., Puckett A. W., Bizyaev 
D., Dilday B., Frieman J. A., Oravetz D. J., Pan K., Quinn T., 
Schneider D. P., Watters S., 2009, ApJ, 695, 268 
Kaib N. A., Quinn T., 2009, Science, 325, 1234 
Krolikowska M., 2001, A&A, 376, 316 
Krolikowska M., 2004, A&A, 427, 1117 
Krolikowska M., 2006, Acta Astronomica, 56, 385 
Krolikowska M., Dybczynski P. A., 2010, MNRAS, 404, 1886 
Levison H. E, Dones L., Duncan M. J., 2001, AJ, 121, 2253 
Marsden B. G., Sekanina Z., Everhart E., 1978, AJ, 83, 64 
Marsden B. G., Sekanina Z., Yeomans D. K, 1973, AJ, 78, 21 1 
Marsden B. G., Williams G. V, 2008, Catalogue of Cometary Or- 
bits 17th Edition. Smithsonian Astrophysical Observatory, Cam- 
bridge, Mass. 
Matese J. J., Lissauer J. J., 2004, Icarus, 170, 508 
Matese J. J., Whitman P. G, 1989, Icarus, 82, 389 
Matese J. J., Whitman P. G., 1992, Celestial Mechanics and Dy- 
namical Astronomy, 54, 13 
Matese J. J., Whitmire D. P., 2011, Icarus, 211, 926 
Morbidelli A., 2005, arXiv:astro-ph/05 12256 
Oort J. H., 1950, Bull.Astron.Inst.Nether., 11, 91 
Oort J. H., Schmidt M., 1951, Bull.Astron.Inst.Nether., 1 1, 259 
Sitarski G, 1989, Acta Astronomica, 39, 345 
Sitarski G, 1998, Acta Astronomica, 48, 547 
Sitarski G, 2002, Acta Astronomica, 52, 471 
Weissman P. R., 1985, in IAU Colloq. 83: Dynamics of Comets: 
Their Origin and Evolution Dynamical Evolution of the Oort 
Cloud, pp 87-96 
Wiegert P., Tremaine S., 1999, Icarus, 137, 84 
Yabushita S„ 1989, AJ, 97, 262 



